A universal Hamiltonian for a quantum dot in the presence of spin-orbit interaction 
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We derive a universal Hamiltonian for a quantum dot in the presence of spin-orbit interaction in 
the symmetry limits of a strong spin-dependent Aharonov-Bohm-like term. We also derive a closed 
expression for the conductance through such a dot, and use it to study the effects of spin-orbit on 
the conductance peak statistics in the presence of an exchange interaction. For a realistic strength of 
the exchange interaction, we find that the width of the peak-spacing distribution is sensitive to spin- 
orbit coupling only in the absence of an orbital magnetic field. We also find that spin-orbit coupling 
modifies the shape of the peak-spacing distribution and suppresses the peak-height fluctuations. 
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The statistical fluctuations of the single-particle spec- 
trum and wavefunctions in a chaotic or diffusive quan- 
tum dot with a large Thouless conductance ot can be 
described by random matrix theory (RMT) 0, 0- In 
open dots, which are strongly coupled to leads, RMT can 
successfully describe the mesoscopic fluctuations of the 
conductance assuming non-interacting electrons. How- 
ever, in almost-isolated dots, which are weakly coupled 
to leads, electron-electron interactions cannot be ignored. 
In such dots, most of the residual interaction terms are 
suppressed in the limit of a large gr, except for a few 
terms. These terms constitute the interacting part of the 
universal Hamiltonian 0-0. and include, in addition to 
the charging energy, a constant exchange interaction. Us- 
ing RMT to describe the single-particle Hamiltonian, a 
significantly better agreement with the data of Refs.lala 
is found once this exchange term is included 0, • In 
general, the exchange interaction suppresses fluctuations 
of both the conductance peak spacings and peak heights. 

Spin-orbit (SO) scattering is also expected to affect 
the conductance fluctuations, although its effects are sup- 
pressed in small dots. Enhanced SO coupling in a parallel 
magnetic field explained the suppression of conductance 
fluctuations in chaotic open 2D GaAs dots 0, [l(J . The 
possible symmetries of the single-particle Hamiltonian in 
such dots with SO scattering were recently classified .11], 
and their signatures were observed in open dots [l2^ . 

Realistic studies of SO effects in almost-isolated dots 
require the inclusion of interactions. However, the uni- 
versal Hamiltonian for a 2D dot has been derived only 
in the absence of SO coupling. Here we study how the 
universal Hamiltonian is modified in the presence of SO 
scattering, both in the presence and absence of an or- 
bital magnetic field. In particular, we derive the universal 
Hamiltonian [Eq. (J5J] in the new symmetry limits intro- 
duced by the leading order SO term (a spin-dependent 
Aharonov-Bohm-like term). We also describe the dot's 
Hamiltonian [Eqs. ©,(171)] in the crossover induced by 
this SO term. Electron correlations in the presence of 
SO scattering were recently studied in metal nanoparti- 
cles |l3j, but such 3D nanoparticles do not possess the 



new symmetries considered here. In addition, we de- 
rive a closed formula for the conductance through a dot 
described by the new universal Hamiltonian [Eqs. © 
- Hll|)]. and study the corresponding conductance peak 
statistics. For a realistic strength of the exchange inter- 
action, we find that the standard deviation of the peak- 
spacing fluctuations is sensitive to SO scattering only in 
the absence of an orbital magnetic field. SO coupling also 
modifies the shape of the peak-spacing distribution and 
suppresses the peak height fluctuations. 

The symmetries of the one-body Hamiltonian in the 
presence of SO in GaAs dots were classified by apply- 
ing a suitable unitary transformation to its original form 
|ll| . The transformed Hamiltonian is expanded in the 
parameter L/X (L is the linear size of the dot and A 
is a mean SO scattering length), which is assumed to 
be small. The two leading contributions are described 
by effective vector potentials that modify the vector po- 
tential of the orbital magnetic field B. The first is a 
spin-dependent Aharonov-Bohm-like vector potential 
of order (L/X) 2 , and the second is a spin- flip term an of 
order (L/X) 3 . In small dots, only the first term is relevant 
and the total effective orbital field is given by 



B c tt — B + B s , 



B so = 4 - (!) 
eX z 



where s z is the electron spin component perpendicular 
to the plane of the dot, and B so is an effective SO field. 
The effect of the SO interaction is described by a dimcn- 
sionless parameter x\ — ngr^so/^o) 2 — k 9t(A/X 2 ) 2 , 
where A is the area of the dot, k is a geometrical coef- 
ficient, $ so is the flux associated with the SO field, and 
$0 = ch/e is the unit flux. In the following we derive the 
universal Hamiltonian in the symmetry limits xj_ 1 
for both B = and B^O. (in practice, the crossover is 
often achieved for x± ~ 1). 

The single-particle eigenstates in the presence of an 
effective field Q are given by \aa), where a = ± de- 
scribes spin up/down electrons with orbital wavefunc- 
tions Va±( r ) an d energies e a ±. In the absence of SO 
coupling (x± — 0), ip a + — tpa—- When SO is present, we 
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have to distinguish between two cases. For B ^ 0, the 
wavefunctions ipa±( r ) correspond to two different values 
B±B BO /2 of the effective field, and become uncorrelated 
for x± 3> 1. Thus, we have a crossover from two de- 
generate Gaussian unitary ensembles (GUE) at x± = 
to two uncorrelated GUE at x± ^> 1 14] . However, for 
-6 = 0, time reversal invariance leads to Kramers degen- 
eracy e a + = e Q _ and ip a -(r) — ipa+{ r )i an d the crossover 
is from two degenerate Gaussian orthogonal ensembles 
(GOE) at x± = to two degenerate GUE at x±_ ^> 1. 

Since \aa) is a complete single-particle basis, we can 
write the dot's Hamiltonian in this basis. We first discuss 
the limits x± = and x± ^> 1 (more precisely, x\ > gx) 
but not the crossover itself. In these limits, we consider 
the "diagonal" part of the interaction, i.e., direct terms 
Vacfycr' \ol(7~^o' j exchange terms ^aa^o'^cran' 1 and Cooper 
channel terms w Q + Q _ ;7 + 7 _. For gx S> 1, we separate this 
"diagonal" interaction into average and fluctuating parts, 
and identify terms that remain finite for gx — * oo. We 



define the following average matrix elements (for 

V$ — V a(J j — a;^cra— a ; ^4 — ^a+a — ;7+7 — ■ (^) 

While all the direct matrix elements have the same av- 
erage V\, we have distinguished two types of exchange 
matrix elements. For x± — 0, the orbital wavefunctions 
are spin-independent and v 2 = v% = J s . For x± —> oo, 
t>2 = J s remains unchanged (the orbital wavefunctions 
correspond to the same spin er), but 1)3 = 0. The van- 
ishing of U3 can be shown, e.g., in a contact model for 
the screened interaction [v(r — r') oc S(r — r')] for which 
oc / c?r^* CT (r)^*_ IT (r)V' 7( T(r)^ Q _ CT (r). For 
B^O, the correlator of the wavefunctions t/j aa and ipa-a 
is a GUE parametric correlator which decays as a power 
law for xi>l F° r 5 = 0, the orbital wavefunc- 

tions have GUE symmetry with ipa-a = il)*L„ and thus 
v 3 oc J dvr aa WVV ( r )</V (r) (r) =0 0. 



I 

Using Eqs. (0), we can write the average part of the diagonal interaction in the form 

Vdiag = ^ {^Vxflafl^ - i« 2 ^ h aa h 1<T + ^V 3 ^ a a< r a i-< r a a-<r<h<' + V 4^T 7 ) - (v r + V 3 + V A ) ^ ^a^a ' ( 3 ) 
a^ 7 a a a 



where h aa is the occupation operator of the state \aa) 
and h a = n a+ + Also TJJ = a' t a+ al l _ is a pair 

creation operator, and we have used w Q ,+ Q ,_ ;Q ,+ Q _ = vi + 
«3 + i>4 (as is easily verified for a contact interaction) . 
The total spin operator of the dot can be represented 

in the basis \aa) as S = ^J2 aT ,aa'( a<7 h CT ') a i<7 cr ^' a "r^ 
(<r are Pauli matrices), and is no longer diagonal in the 
orbital label a. However, since (aa\ja) = <5 Q7 , the z 
component of the total spin is diagonal S z — ^ — 
h a -) = J2 a ^az, and we have ^"mV = k™<**h + 
2s az s lz . The V3 term in Eq. J3J is absent in the limit 
x±_ — > 00, while for x± = it can be related to the spin 
operators through a ! a+ at ( _a a _a~, /+ + a/ a _aL + a a+ a^_ — 
4(s az s 7 2 — s a -s 7 ). Thus in both symmetry limits x±_ = 
and x± — > 00 we can write in the form 

Vdiag = -(wi - v 2 /2)h 2 - i(wi - v 2 )n ~ W3S 2 

- (^-«3)^+«4ftf,(4) 

where T' — J2 a ■ I n the following we denote (3—1 
{(3 = 2) for B = (B ^ 0). For i ± = 0, u 2 = w 3 = Js 
and W4 = 5/31 Jc (Jc is the strength of the Cooper chan- 
nel interaction), and we recover the universal Hamilto- 
nian 0, 0. However, for x± 1, we have v 3 = 
and the S 2 interaction is replaced by S 2 Z . Furthermore, 
«4 = 5p% J c still holds for x±_ ^> 1. For B ^ 0, V4 is obvi- 
ously zero, while for _B = we use = to obtain 

Wq+q-;7+7- K / ^ r |V'a+( r )| 2 |' ! /'7+( r )| 2 an d thus U4 7^ 0. 



As in the absence of SO, the off-diagonal elements of 
the residual interaction are suppressed at large gr- Thus, 
for x± ^> 1, we obtain a new universal Hamiltonian 

H = J2 ^ a n aa + hj d n 2 - J S S 2 Z + SpiJcFf . (5) 

For B = 0, the spin up/spin down levels are degener- 
ate GUE levels (e a + — e a -), while for B ^ 0, the spin 
up/spin down levels are uncorrelated GUE levels. The 
important new feature of JSJ is that the exchange inter- 
action is now given by —J S S Z instead of — J S S 2 [is| . 

In the crossover itself (finite x± ) , the fluctuations of the 
off-diagonal matrix elements are enhanced and cannot 
be ignored Instead, we use the interaction of the 

universal Hamiltonian in the absence of SO scattering, 
and add SO coupling to the single-particle Hamiltonian 

H c = Y,^n aa + hj d n 2 - J S S 2 + SfiiJcftf 

- ^J S (S+S- + S-S+) , (6) 

where we have used S 2 = (S+S- +S-S+)/2 + S 2 . The 
spin and pairing operators can be rewritten in the SO 
eigenstates \aa). The spin projection S z — h X)a — 

n Q _) and the pair operator — J2a a a+ a a- (f° r P = 
1) remain diagonal as in the universal Hamiltonian JSJl. 



3 



However, the spin components in the plane of the dot 
acquire an off-diagonal form 

S+ = ^ C*a a a+ a "r- ! §- = 5ZCa 7 al-a 7+ . (7) 

ay ay 

Here Cay are fluctuating quantities whose statistical prop- 
erties depend on x±. For B^O, ( ay = (a(0)|7(a:j_)) de- 
scribe parametric overlaps in the Gaussian unitary pro- 
cess (GUP) between the eigenstates at scaled parameter 
values x± = and x± 0, 0. For B = 0, Cay are 
just the orthogonal invariants p ai defined by the "real" 
scalar product of the eigenstates a and 7 |l9j. In the 
GUP, K^TF = zi/[(e a " e 7 ) 2 /A 2 + (^i) 2 ], where e, 
are the energy levels at a parameter value x± 3> 1 |lq 
(a s imilar expression holds in the GOE to GUE crossover 
jlflj). Thus the number of levels 7 coupled to a level a 
in (JJJ is ~ Xj_. Since J5J describes an effective Hamilto- 
nian of ~ gx levels around the Fermi energy, it is valid 
for x\ <C <?t- The matrix elements CayC^v of the in- 
teraction S+S- + S-S+ have rms values < l/^xi ano - 
dominate the corrections to JBJl which are of the order 
1/<7t ■ For a;^ > gr, the effective Hamiltonian is the uni- 
versal Hamiltonian The spin S is no longer a good 
quantum number of the Hamiltonians JHJ and Jfyl (for 
^ 0), but S z — M remains a good quantum number. 

In the following we focus on the universal Hamilto- 
nian © and ignore the Cooper channel term. Since 
[n\ a ,S z ] = 0, the occupations n = {ri\ a } form a com- 
plete set of good quantum numbers. The corresponding 
eigenstates \NnM) have energies £^ = ^2\ a £ \a n Xa + 
UdN 2 /2~J S M 2 , where N — n Xa is the total number 
of electrons and M — J2\( n x+ ~ n A-)/2. 

We have calculated the conductance G at temperature 
T using the rate equations approach |20|. Defining a 
scaled conductance g — (KkT / e 2 T)G (f is an average 
tunneling width), we find 

9 = ^wx a gx,y . (8) 

Act 

Here g Xa = (2/f W\ a r x J(T\ a + T\ a ) are the single- 
particle level conductances, where r^ CT (r^) is the par- 
tial width of an electron in level Xa to decay to the left 
(right) lead. The thermal weights w\„ are given by 

w a<t = £ P^m M ) , (9) 

where the sum over all occupation sequences is restricted 
to the level A with spin a being empty. In Eq. |JSJ, f(x) — 
(l + e x / kT )~ 1 is the Fermi-Dirac function evaluated at an 
energy e XM — e\±=fJ s (M ±1/4) —ip, which corresponds 
to the addition of an electron with spin a is to a level A 
(?f is an effective Fermi energy). The quantity P^m = 
e~[ e nM- e F N ]/ kT jZ is the equilibrium probability of the 



state |JVnM), with a partition function Z defined as a 
Boltzmann- weighted sum over all possible N- and (N + 
l)-body states. The sum over all occupation numbers 
in Eq. © can be evaluated in closed form at constant 
M. The constraint n\ a — is taken into account by 
introducing the factor 1 — n\ a , and we have 

w Xa = £(1 - (nxKjPNMfcU) ; ( 10 ) 

M 

where -N/2 < M < N/2. Here n± = N/2 ± M is the 
number of spin- up and spin-down electrons, and Pn,m 
is the probability to find the dot with N electrons and 
spin projection S z = M . Since the trace at fixed N, M 
is equivalent to a trace at fixed n + , n_ , we find 

6 -(F++F- +U N}M )/kT ,„ , . 

Pn,m = e ™+ /Z , (11) 

with U n ,m = U d N 2 /2 - J S M 2 - e F N. In Eqs. 
and (HU F° a = -kT\ntv n<j e-^ >■ e ^ c U>-/ kT and (ri\)° a 
are free energy and canonical occupations of n CT non- 
interacting spinless fermions with energies t\ a . 

Effects of the SO interaction on the conductance peak 
height statistics in the absence of exchange interaction 
were discussed in Ref. 0. Here we consider SO effects 
on both the peak spacing and peak height statistics in 
the presence of exchange. We calculated these statis- 
tics in the limit x± 3> 1 using the universal Hamiltonian 
and compared them with the corresponding statis- 
tics in the absence of SO interaction (i.e., x± — 0) |7j. 
The left panel of Fig. shows the width er(A2) of the 
peak spacing distribution versus kT/A for several val- 
ues J s of the exchange interaction and for both limits 
x± 3> 1 (dashed lines) and x± = (solid lines). In the 
absence of exchange (J s = 0), the SO interaction leads 
to a strong suppression of the spacing fluctuations. How- 
ever, at J s — 0.3A the width is no longer sensitive to the 
SO coupling, except for a small suppression at higher 
temperatures. The peak-spacing distribution itself is af- 
fected by SO coupling and decreases more slowly at large 
spacings. In particular, the bimodality of the x± = dis- 
tribution for kT < 0.3 A disappears for x± ^> 1 (see inset 
in Fig. ^) [2^ • The right panel of Fig. ^ shows similar 
results as in the left panel but for B = 0. In contrast 
to the B ^ case, we observe that the width cr(A 2 ) is 
still sensitive to the SO coupling at J s = 0.3 A. SO en- 
hances the spacing fluctuations at low temperatures but 
suppresses them at higher temperatures. 

A quantity that characterizes the peak height statistics 
is the ratio between the standard deviation cr(g max ) and 
the average <? max of the peak heights (? max . The left panel 
of Fig. |3 shows this ratio as a function of kT/A for J s = 
0.3 A and B ^ in both limits x± = and x± > 1. Wc 
observe that SO scattering suppresses cr(<7 ma x)/ffmax, and 
this suppression becomes stronger with temperature. 

At higher temperatures, inelastic scattering becomes 
important. The right panel of Fig. [21 shows similar re- 
sults as in the left panel but in the rapid-thermalization 
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FIG. 1: The width <r(A2) of the peak-spacing distribution (in 
units of A) for several values J a of the exchange interaction 
and in the presence (left panel) and absence (right panel) 
of an orbital magnetic field. The results in the presence of 
SO coupling for the universal Hamiltonian (dashed lines), 
are compared with the results in the absence of SO coupling 
(solid lines). Inset: the peak-spacing distribution P(A2) (A2 
is measured relative to e 2 /C) in the limits x± = (solid) and 
i 1 »l (dashed) for J s = 0.3A and kT = 0.2A. 



limit of strong inelastic scattering. We observe that 
SO interaction can lead to a significant suppression of 
c((?max)/9max a l s0 m the rapid-thermalization limit. 

At present, no statistical data are available for observ- 
ing SO effects in the presence of an exchange interac- 
tion. Using the value of A determined from experiments 
in large open dots ^3 1 we estimate 11 C 1 for the small 
dots of Refs. 0, 0- It would be interesting to measure the 
conductance peak statistics in almost-isolated dots with 
large area A, in which SO effects are enhanced (for a 
fixed electron density, x± oc A 5 ^ 4 ). In such large dots, it 
is difficult to reach the limit T« A, but our results are 
not restricted to low temperatures. 

In conclusion, we have derived the universal Hamil- 
tonian of a quantum dot in the new symmetry limits 
when the leading order SO interaction term is included, 
both in the presence and absence of an orbital magnetic 
field. Using this universal Hamiltonian, we have identi- 
fied the signatures of SO scattering in the conductance 
peak statistics in the presence of an exchange interac- 
tion. The universal Hamiltonians for other symmetries 
in the presence of SO scattering and the corresponding 
crossover Hamiltonians will be discussed elsewhere. 

We acknowledge useful discussions with B.L. Altshuler, 
P.W. Brouwer, V. Fal'ko, D. Huertas-Hernando, CM. 
Marcus, and A.D. Mirlin. This work was supported in 
part by the U.S. DOE grant No. DE-FG-0291-ER-40608. 
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FIG. 2: The ratio o~(gi max ) /(; ma x between the standard devi- 
ation and the average value of the peak height versus tem- 
perature kT for an exchange interaction of J s = 0.3 A in the 
presence of an orbital field and for both the elastic (left panel) 
and rapid-thermalization (right panel) limits. The dashed 
lines correspond to the universal Hamiltonian while the 
solid lines are the results in the absence of SO coupling 7]. 



[2] 

[3] 

[4] 

[5] 
[6] 
[7] 

[§ 

[9 
[10 

[11 

[12 
[13 

[14 

[15; 

[16 
[17 



[18. 

[19 
[20 

[21 

[22 



[1] Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 



T. Guhr, A. Miiller-Groeling, and H. A. Weidenmuller, 
Phys. Rep. 299, 190 (1998). 

I. L. Kurland, I. L. Aleiner, and B. L. Altshuler, Phys. 
Rev. B 62, 14886 (2000). 

I.L. Aleiner, P.W. Brouwer, and L.I. Glazman, Phys. 
Rep. 358, 309 (2002). 

S. R. Patel et al., Phys. Rev. Lett. 80, 4522 (1998). 
S. R. Patel et al.. Phys. Rev. Lett. 81, 5900 (1998). 
Y. Alhassid and T. Rupp, Phys. Rev. Lett. 91, 056801 
(2003). 

G. Usaj and H. U. Baranger, Phys. Rev. B 67, 121308 
(2003). 

J. A. Folk et al, Phys. Rev. Lett. 86, 2102 (2001). 
B. I. Halperin et al, Phys. Rev. Lett. 86, 2106 (2001). 
I.L. Aleiner and V.I. Fal'ko, Phys. Rev. Lett. 87, 256801 
(2001); 89, 079902(E) (2002). 

D. M. Zumbiihl et al, Phys. Rev. Lett. 89, 276803 (2002). 
D. A. Gorokhov and P.W. Brouwer, Phys. Rev. Lett. 91, 
186602 (2003) ; |cond-mat/031 1086] 

We assume the orbital field to be sufficiently large so that 
none of the effective fields in Q vanish as B so increases. 
Y. Alhassid and H. Attias, Phys. Rev. Lett. 74, 4635 
(1995); Phys. Rev. Lett. 76, 1711 (1996). 
M. Wilkinson and P. Walker, J. Phys. A 28, 6143 (1995). 
The quantity v a a-,-a;~t-aaa is not gauge invariant and 
its average is generally not well defined. However in the 
limit x± — > 00, this quan tity self averages to zero. In 
particular, ^ct^-a^-aa^ 2 — (up to 1/gr corrections). 
The exchange interaction in JSJ is formally similar to the 
dipolar interaction of the double dot model by S. Adam, 
P. W. Brouwer, and P. Sharma, cond-mat/0309074 
S. Adam et al., Phys. Rev. B 66, 165310 (2002). 
Y. Alhassid, T. Rupp, A. Kaminski, and L. I. Glazman, 
arXiv:cond-mat/0212072 

K. Held, E. Eisenberg, and B. L. Altshuler, Phys. Rev. 
Lett. 90, 106802 (2003). 

Finite- gr corrections to the universal Hamiltonian can 
also suppress the bimodality. 



